ChemRAP uncovers specific mRNA translation regulation via RNA 5' phospho-methylation.

5'-end modifications play key roles in determining RNA fates. Phospho-methylation is a noncanonical cap occurring on either 5'-PPP or 5'-P ends. We used ChemRAP, in which affinity purification of cellular proteins with chemically synthesized modified RNAs is coupled to quantitative proteomics, to identify 5'-Pme "readers". We show that 5'-Pme is directly recognized by EPRS, the central subunit of the multisynthetase complex (MSC), through its linker domain, which has previously been involved in key noncanonical EPRS and MSC functions. We further determine that the 5'-Pme writer BCDIN3D regulates the binding of EPRS to specific mRNAs, either at coding regions rich in MSC codons, or around start codons. In the case of LRPPRC (leucine-rich pentatricopeptide repeat containing), a nuclear-encoded mitochondrial protein associated with the French Canadian Leigh syndrome, BCDIN3D deficiency abolishes binding of EPRS around its mRNA start codon, increases its translation but ultimately results in LRPPRC mislocalization. Overall, our results suggest that BCDIN3D may regulate the translation of specific mRNA via RNA-5'-Pme.

Transcriptomics and chromatin accessibility in multiple African population samples.

Mapping the functional human genome and impact of genetic variants is often limited to European-descendent population samples. To aid in overcoming this limitation, we measured gene expression using RNA sequencing in lymphoblastoid cell lines (LCLs) from 599 individuals from six African populations to identify novel transcripts including those not represented in the hg38 reference genome. We used whole genomes from the 1000 Genomes Project and 164 Maasai individuals to identify 8,881 expression and 6,949 splicing quantitative trait loci (eQTLs/sQTLs), and 2,611 structural variants associated with gene expression (SV-eQTLs). We further profiled chromatin accessibility using ATAC-Seq in a subset of 100 representative individuals, to identity chromatin accessibility quantitative trait loci (caQTLs) and allele-specific chromatin accessibility, and provide predictions for the functional effect of 78.9 million variants on chromatin accessibility. Using this map of eQTLs and caQTLs we fine-mapped GWAS signals for a range of complex diseases. Combined, this work expands global functional genomic data to identify novel transcripts, functional elements and variants, understand population genetic history of molecular quantitative trait loci, and further resolve the genetic basis of multiple human traits and disease.

Integrative analyses highlight functional regulatory variants associated with neuropsychiatric diseases.

Noncoding variants of presumed regulatory function contribute to the heritability of neuropsychiatric disease. A total of 2,221 noncoding variants connected to risk for ten neuropsychiatric disorders, including autism spectrum disorder, attention deficit hyperactivity disorder, bipolar disorder, borderline personality disorder, major depression, generalized anxiety disorder, panic disorder, post-traumatic stress disorder, obsessive-compulsive disorder and schizophrenia, were studied in developing human neural cells. Integrating epigenomic and transcriptomic data with massively parallel reporter assays identified differentially-active single-nucleotide variants (daSNVs) in specific neural cell types. Expression-gene mapping, network analyses and chromatin looping nominated candidate disease-relevant target genes modulated by these daSNVs. Follow-up integration of daSNV gene editing with clinical cohort analyses suggested that magnesium transport dysfunction may increase neuropsychiatric disease risk and indicated that common genetic pathomechanisms may mediate specific symptoms that are shared across multiple neuropsychiatric diseases.

Deep learning-assisted genome-wide characterization of massively parallel reporter assays.

Massively parallel reporter assay (MPRA) is a high-throughput method that enables the study of the regulatory activities of tens of thousands of DNA oligonucleotides in a single experiment. While MPRA experiments have grown in popularity, their small sample sizes compared to the scale of the human genome limits our understanding of the regulatory effects they detect. To address this, we develop a deep learning model, MpraNet, to distinguish potential MPRA targets from the background genome. This model achieves high discriminative performance (AUROC = 0.85) at differentiating MPRA positives from a set of control variants that mimic the background genome when applied to the lymphoblastoid cell line. We observe that existing functional scores represent very distinct functional effects, and most of them fail to characterize the regulatory effect that MPRA detects. Using MpraNet, we predict potential MPRA functional variants across the genome and identify the distributions of MPRA effect relative to other characteristics of genetic variation, including allele frequency, alternative functional annotations specified by FAVOR, and phenome-wide associations. We also observed that the predicted MPRA positives are not uniformly distributed across the genome; instead, they are clumped together in active regions comprising 9.95% of the genome and inactive regions comprising 89.07% of the genome. Furthermore, we propose our model as a screen to filter MPRA experiment candidates at genome-wide scale, enabling future experiments to be more cost-efficient by increasing precision relative to that observed from previous MPRAs.

Multiple causal variants underlie genetic associations in humans.

Associations between genetic variation and traits are often in noncoding regions with strong linkage disequilibrium (LD), where a single causal variant is assumed to underlie the association. We applied a massively parallel reporter assay (MPRA) to functionally evaluate genetic variants in high, local LD for independent cis-expression quantitative trait loci (eQTL). We found that 17.7% of eQTLs exhibit more than one major allelic effect in tight LD. The detected regulatory variants were highly and specifically enriched for activating chromatin structures and allelic transcription factor binding. Integration of MPRA profiles with eQTL/complex trait colocalizations across 114 human traits and diseases identified causal variant sets demonstrating how genetic association signals can manifest through multiple, tightly linked causal variants.

BCDIN3D RNA methyltransferase stimulates Aldolase C expression and glycolysis through let-7 microRNA in breast cancer cells.

Type II diabetes (T2D) and specific cancers share many risk factors, however, the molecular mechanisms underlying these connections are often not well-understood. BCDIN3D is an RNA modifying enzyme that methylates specific precursor microRNAs and tRNAHis. In addition to breast cancer, BCDIN3D may also be linked to metabolism, as its gene locus is associated with obesity and T2D. In order to uncover metabolic pathways regulated by BCDIN3D in cancer, we performed an unbiased analysis of the metabolome, transcriptome, and proteome of breast cancer cells depleted for BCDIN3D. Intersection of these analyses showed that BCDIN3D-depleted cells have increased levels of Fructose 1,6 Bisphosphate (F1,6-BP), the last six-carbon glycolytic intermediate accompanied by reduced glycolytic capacity. We further show that elevated F1,6-BP is due to downregulation of Aldolase C (ALDOC), an enzyme that cleaves F1,6-BP mainly in the brain, but whose high expression/amplification is associated with poor prognosis in breast cancer. BCDIN3D regulates ALDOC through a non-canonical mechanism involving the crucial let-7 microRNA family and its target site on the 3'UTR of ALDOC. Overall, our results reveal an important connection between BCDIN3D, let-7 and glycolysis that may be relevant to breast cancer, obesity, and T2D.

Post-transcriptional regulation of antiviral gene expression by N6-methyladenosine.

Type I interferons (IFNs) induce hundreds of IFN-stimulated genes (ISGs) in response to viral infection. Induction of these ISGs must be regulated for an efficient and controlled antiviral response, but post-transcriptional controls of these genes have not been well defined. Here, we identify a role for the RNA base modification N6-methyladenosine (m6A) in the regulation of ISGs. Using ribosome profiling and quantitative mass spectrometry, coupled with m6A-immunoprecipitation and sequencing, we identify a subset of ISGs, including IFITM1, whose translation is enhanced by m6A and the m6A methyltransferase proteins METTL3 and METTL14. We further determine that the m6A reader YTHDF1 increases the expression of IFITM1 in an m6A-binding-dependent manner. Importantly, we find that the m6A methyltransferase complex promotes the antiviral activity of type I IFN. Thus, these studies identify m6A as having a role in post-transcriptional control of ISG translation during the type I IFN response for antiviral restriction.

Evaluating the Genomic Parameters Governing rAAV-Mediated Homologous Recombination.

Recombinant adeno-associated virus (rAAV) vectors have the unique ability to promote targeted integration of transgenes via homologous recombination at specified genomic sites, reaching frequencies of 0.1%-1%. We studied genomic parameters that influence targeting efficiencies on a large scale. To do this, we generated more than 1,000 engineered, doxycycline-inducible target sites in the human HAP1 cell line and infected this polyclonal population with a library of AAV-DJ targeting vectors, with each carrying a unique barcode. The heterogeneity of barcode integration at each target site provided an assessment of targeting efficiency at that locus. We compared targeting efficiency with and without target site transcription for identical chromosomal positions. Targeting efficiency was enhanced by target site transcription, while chromatin accessibility was associated with an increased likelihood of targeting. ChromHMM chromatin states characterizing transcription and enhancers in wild-type K562 cells were also associated with increased AAV-HR efficiency with and without target site transcription, respectively. Furthermore, the amenability of a site to targeting was influenced by the endogenous transcriptional level of intersecting genes. These results define important parameters that may not only assist in designing optimal targeting vectors for genome editing, but also provide new insights into the mechanism of AAV-mediated homologous recombination.

Impact of admixture and ancestry on eQTL analysis and GWAS colocalization in GTEx.

BACKGROUND: Population structure among study subjects may confound genetic association studies, and lack of proper correction can lead to spurious findings. The Genotype-Tissue Expression (GTEx) project largely contains individuals of European ancestry, but the v8 release also includes up to 15% of individuals of non-European ancestry. Assessing ancestry-based adjustments in GTEx improves portability of this research across populations and further characterizes the impact of population structure on GWAS colocalization. RESULTS: Here, we identify a subset of 117 individuals in GTEx (v8) with a high degree of population admixture and estimate genome-wide local ancestry. We perform genome-wide cis-eQTL mapping using admixed samples in seven tissues, adjusted by either global or local ancestry. Consistent with previous work, we observe improved power with local ancestry adjustment. At loci where the two adjustments produce different lead variants, we observe 31 loci (0.02%) where a significant colocalization is called only with one eQTL ancestry adjustment method. Notably, both adjustments produce similar numbers of significant colocalizations within each of two different colocalization methods, COLOC and FINEMAP. Finally, we identify a small subset of eQTL-associated variants highly correlated with local ancestry, providing a resource to enhance functional follow-up. CONCLUSIONS: We provide a local ancestry map for admixed individuals in the GTEx v8 release and describe the impact of ancestry and admixture on gene expression, eQTLs, and GWAS colocalization. While the majority of the results are concordant between local and global ancestry-based adjustments, we identify distinct advantages and disadvantages to each approach.

Transcriptomic signatures across human tissues identify functional rare genetic variation.

Rare genetic variants are abundant across the human genome, and identifying their function and phenotypic impact is a major challenge. Measuring aberrant gene expression has aided in identifying functional, large-effect rare variants (RVs). Here, we expanded detection of genetically driven transcriptome abnormalities by analyzing gene expression, allele-specific expression, and alternative splicing from multitissue RNA-sequencing data, and demonstrate that each signal informs unique classes of RVs. We developed Watershed, a probabilistic model that integrates multiple genomic and transcriptomic signals to predict variant function, validated these predictions in additional cohorts and through experimental assays, and used them to assess RVs in the UK Biobank, the Million Veterans Program, and the Jackson Heart Study. Our results link thousands of RVs to diverse molecular effects and provide evidence to associate RVs affecting the transcriptome with human traits.

Genetic regulation of gene expression and splicing during a 10-year period of human aging.

BACKGROUND: Molecular and cellular changes are intrinsic to aging and age-related diseases. Prior cross-sectional studies have investigated the combined effects of age and genetics on gene expression and alternative splicing; however, there has been no long-term, longitudinal characterization of these molecular changes, especially in older age. RESULTS: We perform RNA sequencing in whole blood from the same individuals at ages 70 and 80 to quantify how gene expression, alternative splicing, and their genetic regulation are altered during this 10-year period of advanced aging at a population and individual level. We observe that individuals are more similar to their own expression profiles later in life than profiles of other individuals their own age. We identify 1291 and 294 genes differentially expressed and alternatively spliced with age, as well as 529 genes with outlying individual trajectories. Further, we observe a strong correlation of genetic effects on expression and splicing between the two ages, with a small subset of tested genes showing a reduction in genetic associations with expression and splicing in older age. CONCLUSIONS: These findings demonstrate that, although the transcriptome and its genetic regulation is mostly stable late in life, a small subset of genes is dynamic and is characterized by a reduction in genetic regulation, most likely due to increasing environmental variance with age.

BCDIN3D regulates tRNAHis 3' fragment processing.

5' ends are important for determining the fate of RNA molecules. BCDIN3D is an RNA phospho-methyltransferase that methylates the 5' monophosphate of specific RNAs. In order to gain new insights into the molecular function of BCDIN3D, we performed an unbiased analysis of its interacting RNAs by Thermostable Group II Intron Reverse Transcriptase coupled to next generation sequencing (TGIRT-seq). Our analyses showed that BCDIN3D interacts with full-length phospho-methylated tRNAHis and miR-4454. Interestingly, we found that miR-4454 is not synthesized from its annotated genomic locus, which is a primer-binding site for an endogenous retrovirus, but rather by Dicer cleavage of mature tRNAHis. Sequence analysis revealed that miR-4454 is identical to the 3' end of tRNAHis. Moreover, we were able to generate this 'miRNA' in vitro through incubation of mature tRNAHis with Dicer. As found previously for several pre-miRNAs, a 5'P-tRNAHis appears to be a better substrate for Dicer cleavage than a phospho-methylated tRNAHis. Moreover, tRNAHis 3'-fragment/'miR-4454' levels increase in cells depleted for BCDIN3D. Altogether, our results show that in addition to microRNAs, BCDIN3D regulates tRNAHis 3'-fragment processing without negatively affecting tRNAHis's canonical function of aminoacylation.

Genetic analyses of human fetal retinal pigment epithelium gene expression suggest ocular disease mechanisms.

The retinal pigment epithelium (RPE) serves vital roles in ocular development and retinal homeostasis but has limited representation in large-scale functional genomics datasets. Understanding how common human genetic variants affect RPE gene expression could elucidate the sources of phenotypic variability in selected monogenic ocular diseases and pinpoint causal genes at genome-wide association study (GWAS) loci. We interrogated the genetics of gene expression of cultured human fetal RPE (fRPE) cells under two metabolic conditions and discovered hundreds of shared or condition-specific expression or splice quantitative trait loci (e/sQTLs). Co-localizations of fRPE e/sQTLs with age-related macular degeneration (AMD) and myopia GWAS data suggest new candidate genes, and mechanisms by which a common RDH5 allele contributes to both increased AMD risk and decreased myopia risk. Our study highlights the unique transcriptomic characteristics of fRPE and provides a resource to connect e/sQTLs in a critical ocular cell type to monogenic and complex eye disorders.

Crosstalk between the RNA Methylation and Histone-Binding Activities of MePCE Regulates P-TEFb Activation on Chromatin.

RNAP II switching from the paused to the productive transcription elongation state is a pivotal regulatory step that requires specific phosphorylations catalyzed by the P-TEFb kinase. Nucleosolic P-TEFb activity is inhibited by its interaction with the ribonuclear protein complex built around the 7SK small nuclear RNA (7SK snRNP). MePCE is the RNA methyltransferase that methylates and stabilizes 7SK in the nucleosol. Here, we report that MePCE also binds chromatin through the histone H4 tail to serve as a P-TEFb activator at specific genes important for cellular identity. Notably, this histone binding abolishes MePCE's RNA methyltransferase activity toward 7SK, which explains why MePCE-bound P-TEFb on chromatin may not be associated with the full 7SK snRNP and is competent for RNAP II activation. Overall, our results suggest that crosstalk between the histone-binding and RNA methylation activities of MePCE regulates P-TEFb activation on chromatin in a 7SK- and Brd4-independent manner.

Chromatin Regulates Genome Targeting with Cisplatin.

Cisplatin derivatives can form various types of DNA lesions (DNA-Pt) and trigger pleiotropic DNA damage responses. Here, we report a strategy to visualize DNA-Pt with high resolution, taking advantage of a novel azide-containing derivative of cisplatin we named APPA, a cellular pre-extraction protocol and the labeling of DNA-Pt by means of click chemistry in cells. Our investigation revealed that pretreating cells with the histone deacetylase (HDAC) inhibitor SAHA led to detectable clusters of DNA-Pt that colocalized with the ubiquitin ligase RAD18 and the replication protein PCNA. Consistent with activation of translesion synthesis (TLS) under these conditions, SAHA and cisplatin cotreatment promoted focal accumulation of the low-fidelity polymerase Polη that also colocalized with PCNA. Remarkably, these cotreatments synergistically triggered mono-ubiquitination of PCNA and apoptosis in a RAD18-dependent manner. Our data provide evidence for a role of chromatin in regulating genome targeting with cisplatin derivatives and associated cellular responses.

Click Quantitative Mass Spectrometry Identifies PIWIL3 as a Mechanistic Target of RNA Interference Activator Enoxacin in Cancer Cells.

Enoxacin is a small molecule that stimulates RNA interference (RNAi) and acts as a growth inhibitor selectively in cancer but not in untransformed cells. Here, we used alkenox, a clickable enoxacin surrogate, coupled with quantitative mass spectrometry, to identify PIWIL3 as a mechanistic target of enoxacin. PIWIL3 is an Argonaute protein of the PIWI subfamily that is mainly expressed in the germline and that mediates RNAi through piRNAs. Our results suggest that cancer cells re-express PIWIL3 to repress RNAi through miRNAs and thus open a new opportunity for cancer-specific targeting.

Small RNA Sequencing in Cells and Exosomes Identifies eQTLs and 14q32 as a Region of Active Export.

Exosomes are small extracellular vesicles that carry heterogeneous cargo, including RNA, between cells. Increasing evidence suggests that exosomes are important mediators of intercellular communication and biomarkers of disease. Despite this, the variability of exosomal RNA between individuals has not been well quantified. To assess this variability, we sequenced the small RNA of cells and exosomes from a 17-member family. Across individuals, we show that selective export of miRNAs occurs not only at the level of specific transcripts, but that a cluster of 74 mature miRNAs on chromosome 14q32 is massively exported in exosomes while mostly absent from cells. We also observe more interindividual variability between exosomal samples than between cellular ones and identify four miRNA expression quantitative trait loci shared between cells and exosomes. Our findings indicate that genomically colocated miRNAs can be exported together and highlight the variability in exosomal miRNA levels between individuals as relevant for exosome use as diagnostics.

MiR-191 Regulates Primary Human Fibroblast Proliferation and Directly Targets Multiple Oncogenes.

miRNAs play a central role in numerous pathologies including multiple cancer types. miR-191 has predominantly been studied as an oncogene, but the role of miR-191 in the proliferation of primary cells is not well characterized, and the miR-191 targetome has not been experimentally profiled. Here we utilized RNA induced silencing complex immunoprecipitations as well as gene expression profiling to construct a genome wide miR-191 target profile. We show that miR-191 represses proliferation in primary human fibroblasts, identify multiple proto-oncogenes as novel miR-191 targets, including CDK9, NOTCH2, and RPS6KA3, and present evidence that miR-191 extensively mediates target expression through coding sequence (CDS) pairing. Our results provide a comprehensive genome wide miR-191 target profile, and demonstrate miR-191's regulation of primary human fibroblast proliferation.

miR-503 represses human cell proliferation and directly targets the oncogene DDHD2 by non-canonical target pairing.

BACKGROUND: The pathways regulating the transition of mammalian cells from quiescence to proliferation are mediated by multiple miRNAs. Despite significant improvements in our understanding of miRNA targeting, the majority of miRNA regulatory networks are still largely unknown and require experimental validation. RESULTS: Here we identified miR-503, miR-103, and miR-494 as negative regulators of proliferation in primary human cells. We experimentally determined their genome wide target profiles using RNA-induced silencing complex (RISC) immunoprecipitations and gene expression profiling. Analysis of the genome wide target profiles revealed evidence of extensive regulation of gene expression through non-canonical target pairing by miR-503. We identified the proto-oncogene DDHD2 as a target of miR-503 that requires pairing outside of the canonical 5' seed region of miR-503, representing a novel mode of miRNA-target pairing. Further bioinformatics analysis implicated miR-503 and DDHD2 in breast cancer tumorigenesis. CONCLUSIONS: Our results provide an extensive genome wide set of targets for miR-503, miR-103, and miR-494, and suggest that miR-503 may act as a tumor suppressor in breast cancer by its direct non-canonical targeting of DDHD2.

A Myc-microRNA network promotes exit from quiescence by suppressing the interferon response and cell-cycle arrest genes.

The transition of mammalian cells from quiescence to proliferation is accompanied by the differential expression of several microRNAs (miRNAs) and transcription factors. However, the interplay between transcription factors and miRNAs in modulating gene regulatory networks involved in human cell proliferation is largely unknown. Here we show that the miRNA miR-22 promotes proliferation in primary human cells, and through a combination of Argonaute-2 immunoprecipitation and reporter assays, we identified multiple novel targets of miR-22, including several cell-cycle arrest genes that mediate the effects of the tumor-suppressor p53. In addition, we found that miR-22 suppresses interferon gene expression by directly targeting high mobility group box-1 and interferon regulatory factor (IRF)-5, preventing activation of IRF3 and NF-κB, which are activators of interferon genes. The expression of interferon genes is elevated in quiescent cells and their expression is inhibitory for cell proliferation. In addition, we find that miR-22 is activated by the transcription factor Myc when quiescent cells enter proliferation and that miR-22 inhibits the Myc transcriptional repressor MXD4, mediating a feed-forward loop to elevate Myc expression levels. Our results implicate miR-22 in downregulating the anti-proliferative p53 and interferon pathways and reveal a new transcription factor-miRNA network that regulates the transition of primary human cells from quiescence to proliferation.